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I review many-body effects on the resistivity of a multiorbital system beyond Landau’s 
Fermi-liquid (FL) theory. Landau’s FL theory succeeds in describing electronic proper¬ 
ties of some correlated electron systems at low temperatures. However, the behaviors 
deviating from the temperature dependence in the FL, non-FL-like behaviors, emerge 
near a magnetic quantum-critical point. These indicate the importance of many-body 
effects beyond Landau’s FL theory. Those effects in multiorbital systems have been little 
understood, although their understanding is important to deduce ubiquitous properties 
of correlated electron systems and characteristic properties of multiorbital systems. To 
improve this situation, I formulate the resistivity of a multiorbital Hubbard model us¬ 
ing the extended Eliashberg theory and adopt this method to the inplane resistivity 
of quasi-two-dimensional paramagnetic ruthenates in combination with the fluctuation- 
exchange approximation including the current vertex corrections arising from the self¬ 
energy and Maki-Thompson term. The results away from and near the antiferromagnetic 
quantum-critical point reproduce the temperature dependence observed in Si^RuCL and 
Sr2Ruo.075Tio.025 04, respectively. I highlight the importance of not only the momentum 
and the temperature dependence of the damping of a quasiparticle but also its orbital 
dependence in discussing the resistivity of correlated electron systems. 

Keywords : many-body effects; non-Fermi-liquid-like behaviors; ruthenates; t2g orbital; 
nearly magnetic metal; fluctuation-exchange approximation; current vertex correction. 


1. Introduction 

Landau’s Fermi-liquid (FL) theorjQ 2] can describe electronic properties of some 
correlated electron systems at low temperatures^ In this theory, low-energy exci¬ 
tations are described by quasiparticles (QPs), i.e. electrons in a self-consistent field 
of surrounding electrons due to electron correlation. Also, the interactions between 
QPs, described by the Landau parameters, are independent of temperature. Due to 
these two properties, the temperature dependences of physical quantities in low-T 
region are governed by the temperature dependence of the Fermi distribution func¬ 
tion, and the corrections due to electron correlation, the mass enhancement and 
the FL correction, are independent of temperature. As a result, thermodynamic or 
magnetic or transport quantities show the same temperature dependences as those 
of the free electron system in low-T region, and the difference is the coefficient due 


1 








June 17, 2015 0:42 WSPC/INSTRUCTION FILE review-rho'final 


2 Naoya Arakawa 


to the mass enhancement or the FL correction or both. For example, the specific 
heat is proportional to T, and the coefficient is renormalized by the mass enhance¬ 
ment; the spin susceptibility is independent of temperature, and the coefficient is 
renormalized by the mass enhancement and the FL correction. 

The original Landau phenomenological theory^ can be justified by microscopic 
perturbation theory with several imposed conditions^. 

One of the basic assumptions, the one-to-one correspondence, becomes valid if 
the QP damping is much smaller than temperature considered (i.e., the QP lifetime, 
the inverse of the QP damping, is very long). Actually, the single-particle spectral 
function near the Fermi level becomes delta-function-type in the coherent limit [i.e., 
la{k)/T -A- 0 for all Fermi momenta], 


A a (k) = --ImG^ R) (fc) = 


,(fc) 


7 a ( fe ) 


7 T [e-£a( fc )] 2 +7aO) 2 


z a (k)S(e-Cm . (1) 


Here we consider only the coherent part, G^(k) = (k)+i-y* (k) ' where the self¬ 

energy is approximated as T^\k) « Ei R ^(fc,0) + e d ^ a 9l , ( ' k) | e ->o; k is k = (fe,e), 
a is the QP band index, £*(fe) is the QP energy, which is of the order of T, 


z a (k) = [1 - 


dReS^ (k,cj) I 


duj 


w _> 0 ] 1 is the inverse of the mass enhancement factor, 


and 7 *(fc) = — z Q (fc)ImSo ’ (k, 0) is the QP damping. Since the delta-function-type 
spectral function is obtained for an exact eigenstate, the FL becomes an approxi¬ 
mate eigenstate for momenta near the Fermi level if the QP dampings at these mo¬ 
menta are much smaller than T. In the FL theory the QP damping becomes small 
rapidly compared with decreasing temperature as a result of its T 2 dependent 0 ®. 

In addition, the other basic assumption about the Landau parameters becomes 
valid if the reducible four-point vertex functions and mass enhancement factor are 
independent of temperaturePEl. Note, first, that the Landau parameter is propor¬ 
tional to the product of the reducible fou r-po int vertex function and the square 
root of the four mass enhancement factors^" second, that the four-point vertex 
function describes the multiple scattering of an electron and a liolePl, third, that 
all reducible diagrams can be split into two parts by removing a pair of the single¬ 
particle Green’s functions, while all irreducible diagrams cannot be done by using 
that removing^. 

In contrast to successful description^ in the paramagnetic (PM) phase of 
Si' 2 Ru 04 , Landau’s FL theory fails in describing electronic properties of other 
PM ruthenates near a magnetic quantum-critical point (QCP). For example, in 
Si' 2 Ruo.o 75 Tio.o 25 G 4 , located near an antiferromagnetic (AF) QCF^, the spin sus¬ 
ceptibility shows the Curie-Weiss-like temperature dependence and the inplane re¬ 
sistivity, pab, shows the T-linear dependence^], deviating from the FL-type T 2 
dependence; the characteristic wave vectors of this AF QCP^, q ss (^, qf) and 
its symmetrically equivalent ones, are the same as the wave vectors of the most 
strongly enhanced spin fluctuation in Sr 2 RuoJEII Also, m C& 2 — x Sr^R/nO^ around 
x = 0.5, located near a ferromagnetic QCP 11 ) the Curie-Weiss-like temperature 
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dependence of the spin susceptibility and T 1A dependence of p a b are observed^. 

Since such non-FL-like behaviors are observed in other systems near a magnetic 
QCp43 01 - a Mott transition^!, where many-body effects generally become very 
important, these experimental factd^EB; indicate the necessity of both considering 
how the basic assumptions of Landau’s FL theory are violated and discussing many- 
body effects beyond Landau’s FL theory in a more elaborate theory. 

So far, there are two candidates for the origin of such non-FL-likc behaviors. 

One is bad metaff^ due to local correlation enhanced near a Mott transition. 
If the low-energy excitations can be satisfactorily described by using only low-order 
Taylor series of the self-energy in terms of frequency, the coherent part of the single¬ 
particle Green’s function plays dominant roles in discussing electronic properties. 
However, if there are some non-negligible contributions from the higher-order Tay¬ 
lor series, the incoherent part also becomes important. Such case is realized in a 
PM metallic phase near the Mott transition due to the formation of the upper 
and lower Hubbard peaks, arising from enhanced local correlation^ Since per¬ 
turbation theory can treat the coherent part appropriately and its treatment of 
the incoherent part is insufficient ED perturbation theory is unsuitable in the lat¬ 
ter case. Instead, the latter case can be well described in dynamical-mean-field 
theory (DMFT) since the DMFT can t ake account of the frequency dependence 
of the self-energy nonperturbatively ED. Actually, several non-FL-like behaviors 
are obtained in the DMFT near the Mott transition as a result of the T-linear 
unrenor mal ized QP damping and the temperature-dependent mass enhancement 
factoJ^SlUJ. [In this paper, I define the unrenormalized QP damping as the QP 
damping divided by the inverse of the mass enhancement factor.] These results in¬ 
dicate the importance of nonpertur bative effects of local correlation near the Mott 
transition. Note that Hund’s metap2123| j s classified as the bad metal since the 
non-FL-like behaviors in the Hund’s metal arise from local spin fluctuation en¬ 
hanced near a Mott transition, although there is a crucial difference between the 
Hund’s rule coupling dependence of the critical value of the infrao rbit al Coulomb 
interaction for the Mott transition at half-filling and non-half-filling^'-’. 

The other is nearly magnetic metaP^ES 26 27 due S p a ti a l correlation en¬ 
hanced near a magnetic QCP. If a system is located far away from a magnetic 
QCP, all scattering processes are independent of temperature. However, if the sys¬ 
tem approaches the QCP, several characteristic scattering processes of the QCP 
show the strong temperature-dependent enhancement^. For example, in case near 
a stripe-type AF QCP the scattering processes mediated by AF spin fluctuations 
with q = (7r, 7 r) and i ts s ymmetrically equivalent ones are strongly enhanced as 
decreasing temperature^. Such strong temperature-dependent enhancement leads 
to the strong temperature dependence of the reducible four-point vertex function 
whose momentum is characterized by the QCP. Thus, the basic assumption of Lan¬ 
dau’s FL theory about the Landau parameters is violated only for the characteristic 
momenta of spin fluctuation enhanced near the QCI^Sl. Also, the basic assump- 
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tion about the QP damping is violated for the characteristic momenta due to the 
formation of hot spot, arising from enhanced spin fluctuation^; at the hot spot, 
the QP damping does not become much smaller than temperature considered. For 
example, in case near the stripe-type AF QCP, the QP dampings at the momenta 
connected by the nesting vector q = (n, ir) and its symmetrically equivalent ones are 
more strongly enhanced than those at the other momenta due to the enhancement 
of the corresponding AF spin fluctuations^ 8 . These two violations suggest the ne¬ 
cessity of discussing electronic properties near a magnetic QCP in the microscopic 
perturbation theory beyond Landau’s FL theory. Actually, fluctuation-exchange 
(FLEX) approximation reproduces several non-FL-like behaviors due to the hot¬ 
spot structure of the QP da mping or the Curie-Weiss-like temperature dependence 
of spin fluctuations or botlJ28 29. a s explained in Sect. 2.3, this approximation 
can take account of spatial correlation beyond a mean-field theory and describe 
electronic properties of a metallic p hase a t low temperatures for moderately strong 
electron correlation satisfactorily 3 ® 3 - 3 21. I n contrast to case near the Mott transi¬ 
tion, the DMFT is inappropriate to describe electronic prope rtie s near a magnetic 
QCP since the DMFT neglects spatial correlation completely321. These results in¬ 
dicate powerfulness of the microscopic perturbation theory and the importance of 
temperature-dependent spatial correlation near a magnetic QCP. 

With the backgrounds explained above, I studied electronic structure and mag¬ 
netic and transport propertied® of ruthenates near and away from the AF QCP 
in the FLEX approximation with current vertex corrections (CVCs) arising from 
the self-energy and Maki-Thompson (MT)^ 3 ^ term for the ^g-orbital Hubbard 
model on a square lattice and succeeded in reproducing several experimental re¬ 
sults of Sr 2 Ruo. 075 Tio .025 OJMEI and Sr 2 RuOi ^ 6 l 37 ^ 8 l. Thus, the non-FL-like 
behaviors 8 ® in Sr 2 Ru 0 . 075 Ti 0 . 02 . 5 O 4 can be understood as the nearly magnetic 
meal near the AF QCP. Moreover, since the results away from the AF QCP about 
the orbital dependence of the mass enhancement are in better agreement with 
the experiment 3 ® in Sr 2 Ru 04 than those of the DMFT 2^1, electronic properties of 
ruthenates at low temperatures except a few cased® 1 ® near the Mott transition 
may be better described in the microscopic perturbation theory than in the DMFT. 

‘ M 


In this paper, I review part of the above previous study ^ 1 and show some 
new results. In Sect. 2, I explain the microscopic theory^! US ed for analyzing p a b 
of some quasi-two-dimensional (quasi-2D) PM ruthenates. In Sect. 3, I show the 
results about many-body effects on p a b of the ruthenates near and away from the 
AF QCP in the FLEX approximation with the CVCs arising from the self-energy 
and MT term and compare these results with the results obtained in other cases 
where the CVCs are more simplified. In Sect. 4, I summarize the results and draw 
some conclusions. 

Due to the limit of space, I do not consider the Aslamasov-Larkin (AL) CVC@2], 
which is the other CVC in the FLEX approximation; its detailed derivation and 
effects are going to be discussed elsewhere 
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2. Method 

In this section, I explain an effective model of some quasi-2D ruthenates, briefly 
review the formal derivation of the resistivity of a multiorbital Hubbard model in 
a PM metallic phase, and formulate the microscopic perturbation theory used to 
calculate the resistivity. The more detailed exp lanations about those derivation and 
formulation are going to be given elsewhere . 

In the following, I use the unit h = c = e = /ib = = 1, set the coordinates x, 

y , and z in the directions of the Ru—O bonds of a RuOe octahedral, and label the 
d xz , d yz , and d xy orbitals as 1, 2, and 3, respectively. 

2.1. Effective model for some quasi-2D ruthenates 

I explain an effective model to describe electronic properties of some quasi-2D 
ruthenates without the rotation or the tilting of RuOg octahedral 

Before introducing the effective model, I briefly explain several basic elec¬ 
tronic properties. Some ruthenates whose crystal structures are 214-type such as 
S^RuOj^or 327-type such as Sr 3 Ru20?^^are categorized into quasi-2D t 2s -orbital 
systems. For simplicity, we focus on Sr 2 Ru 04 ; the following properties remain qual¬ 
itatively the same in other metallic ruthenateJ^ First, the inplane resistivity is 
about 10 -3 times as small as the out-of-plane resistivity at low temperatures^®!, and 
the almost cylindrical Fermi surface (FS) is observed in the de Haas-van Alphen 
measurement®®. These indicate quasi-2D electronic conduction. Moreover, accord¬ 
ing to several density-functional calculations^^ in local-density approximation 
(LDA), conducting bands near the Fermi level are formed by the antibonding bands 
of the Ru t 2 g and the O 2 p orbitals, and the t 2 g orbitals mainly contribute to the 
density-of-states (DOS) near the Fermi level. Since the topology of the FS obtained 
in the LDA qualitatively agrees with experiment the orbitals play domi¬ 
nant roles in discussing electronic properties at low temperatures. 

With the above background, I assume that the electronic structure obtained in 
the LDA™ 

for Sr 2 Ru 04 is a good starting point to consider many-body effects 
beyond a mean-field approximation, and I use a ^g-orbital Hubbard model on a 
square lattice as the effective model. Thus, the Hamiltonian is H = Hq + H m t with 


Ho “E E EE e a,b{k)c\. as Ckbs 


( 2 ) 


k a,b— 1 s=t,i 


and 


Hint = jE E 

^ n n h r rl 


j a,b : c,d—l si ,S 2 ,S3 ,S4="f,4, 
3 


TTSlS 2 S 3 S 4 -t -t * Z 
u abed c jas 1 c jds4 C 3 cs 3 c 3bs 2 


3 


^ EE EE + U JE EE EE ftjaftjb 


j a= 1 


j a—1b<a 


*^H E E E (25ja ‘ Sjb + 2 fi'jofi'jb) + J EE EE EE aX^jbi^jbf’i (3) 


3 
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Fig. 1. (a) Band structure, (b) FS, and (c) DOS for Hq whose parameters are chosen so as to 
reproduce the results in the LDA 47,48 for Sr 2 Ru 04 . 


where e a b(k) is the energy dispersions, measuring from the chemical potential, /i, 
U, U ', Jh, and J 7 are intraorbital Coulomb interaction, interorbital Coulomb in¬ 
teraction, Hund’s rule coupling, and pair hopping term, hj a is hj a = ^2 s fij as = 
E s CjasCjas and is Sj a = \ E s s’ ^jas^^'Cjas 1 with the Pauli matrices <T SS >. 
Since I do not consider the effectJ^Si 0 f the rotation and the tilting, we focus on 
electronic properties of Sr 2 RuO,p- and some doped S^RuO^ without these distor¬ 
tions. Note that the rotation is present in Ca 2 -xSr :E Ru 04 around x = 0.5^. 

By considering some symmetrically possible hopping processes and the difference 
between the crystalline-electric-field energies of the d xy and d xz / yz orbitals, A t 2g , 
we can construct the tight-binding model whose e a b{k) is given by 


en(fc) = - 


A. 


*2g 


— 2 ti cos k x — 2 f 2 cos k y — n, 


ei 2 (fc) = e 2 i(fc) = At' sin k x sin k y , 

e 2 2 (fe) =- 5 ^- — 2 1, 2 cos k x — 2 ti cos k y — /j,, 


£33(fe) = 


2A 


*2g 


— 2f 3 (cos k x + cos k y ) — At 4 cos k x cos k y — jj,, 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


and otherwise e ab (k) = 0 . n is chosen so that the electron number per a site, n e , is 
fixed. For the actual calculations, fi is determined by the bisection method using 


9 2 T 3 r 

ne = -Tf EE/( e «( fe )) + -yEEE G aa (fc,iw m ) - G° aa (k,iu m ) 

^ k a M k m a =1 


( 8 ) 


where the second term becomes zero without the interaction terms. Here e a (k) is 


.(*) = E ( U° k \ a e ab (k)(U ° k ) ba , 

a,6=1 


( 9 ) 


/(e) is the Fermi distribution function, G° b (fc,iu; m ) is the noninteracting single¬ 
particle Green’s function with fermionic Matsubara frequency, u> m = nT(2m +1), 

G° ab (k,ico m ) = E(u° k ) aa -. —^7r(£4 0t U (io) 

and G a b(k,iui m ) is the single-particle Green’s function whose determination is 
explained in Sect. 2.3. In Eq. ©, we put the chemical potentials in f(e a (k)), 
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G aa (fc,iw m ), and G° a (fc,iw m ) the same to reduce the numerical error arising 
from the cut-off frequency. To reproduce the electronic structure obtained in the 
LDA 14ILi8] Iget , At 29 ) = (0.675,0.09,0.45,0.18,0.03,0.13) (eV) and 

n e = 4. Actually, we see from Fig. |T] that the total bandwidth, the topology of 
each FS sheet, and the location of the van Hove singularity of the d xy orbital agree 
with the LDA resultsP^HEl the total bandwidth is about 4 eV; the FS consists of 
the quasi-ID hole-like a and electron-like /3 sheets of the d xz / yz orbital and the 
quasi-2D electron-like 7 sheet of the d xy orbital, and the 7 sheet is located nearer 
the inner sheet in k x = k y or k x = — k y line than in the experiment"^® the van 
Hove singularity is located above the Fermi level. Also, the occupation numbers of 
the d xz / yz and the d xy orbital, being n xz / yz = 1.38 and n xy = 1.25, are consistent 
with the LDA values^ 3HI. 

Then, I set J' = Jr, U' = U — 2Jh, and J H = y-, use U as a parameter, and 
treat the effects of interactions in the FLEX approximation with the CVCs arising 
from the self-energy and MT term; its detail is explained in Sect. 2.3. 

Finally, we remark on suitability neglecting the spin-orbit coupling of Ru ions for 
discussing many-body effects on the resistivity. A density-functional calculation^ 
for S^RuCH within local-spin-density approximation shows that the coupling con¬ 
stant is 0.167 eV, and that the main effect on the electronic structure is the weak 
mixing between the bands of the d xz / yz and the d xy orbital around k = (yy,yy) and 
its symmetrically equivalent ones. This effect will not qualitatively change the re¬ 
sults shown in Sect. 3 since this coupling constant is smaller than the main terms in 
H and that weak mixing will lead to small changes of the momentum dependence 
and value of the spin susceptibility from those without the spin-orbit coupling. 
Thus, I believe that neglecting the spin-orbit coupling is suitable for qualitative 
discussions about many-body effects on the resistivity. 


2.2. Extended Eliashberg theory for the resistivity of a multiorbital 
Hubbard model in a PM metallic phase 

I briefly review the formal derivation"^ ^3‘ 0 f the resistivity of a multiorbital Hub¬ 
bard model in a PM metallic phase. We first derive an exact expression of the 
longitudinal conductivity, the inverse of the resistivity, in the presence of electron 
correlation within the linear-response theory^. Then, we rewrite this exact ex¬ 
pression in terms of the four-point vertex function by using the three-point vector 
vertex function. Due to difficulty solving the exact expression, we derive an ap¬ 
proximate expression by using the most-divergent-term approximation introduced 
by Eliashberg . 

To discuss the resistivity within the linear-response theory, we use the Kubo 
formula^ for the longitudinal conductivity, ov„ (u = x, y), in w-limit and a>Tt ra ns 
1 with t trans being the transport relaxation time^l, which is of the order of magni- 
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tude of the QP lifetime. Namely, is given by 


er„ 2 , = 2 lim lim 

uj —>0 q—t 0 


0 ) 


= 2 lim 

cj—>-0 


K™ 

l-'-l'V 


( 0 , 


UJ 


— 

i-*- vv 


( 0 , 0 ) 


( 11 ) 


where K^ (0,w) is obtained by the analytic continuation of K uu (ifl n ), 

I\„ v (iQ, n ) = lim-i: J Q T dre lQ " T (T r J q ^(r) J_q„(0)} 

= 4 E E/ 0 T dre in - r (nfc y ) 6 aKv)c<i(T T cl 6 (r)c fc a(r)cl'c 2 fc'd) 

-' V k,k' {a} 

= "77 X (‘^ku ) hn. (f-'k 1 u )r/l Kg hr.djk ■ fc (12) 

fc.fc' {a} 


with bosonic Matsubara frequency, f2„ = 2nTn. Here X){a} I s X){a} — X)abcd’ 
(' v kv )ab is the group velocity, ( v kv ) a b = d£ g b k [ k) , and K abcd (k, k'; iQ n ) is 


K-abcdik^h 5 — $k,k'T G ac (k^ ^^m+ri)G db{k i ^m) 

m 

T GaA^fe-t ^m-\-n)GdD{k , ^m) 

m,m' {A} 

X Gcc{k i ^m'-\-n)^{A}{k^ iid mi k , iiOrn' ] Al n ), (U) 


where Tr^Gfe, iu) m , k’, zw m /; *f2„) = r ABCD{k,iujm,k' ,iui rn i]iQ, n ) is the reducible 
four-point vertex function. Thus, the analytic continuation of K a b c d{k , fc'; zH„) is 
necessary to calculate a^ v . 

Before the analytic continuation of K a bcd(k , kiCl n ), I remark on the important 
physical meanings of lim^^o li m q->o and WT tran s "C 1. The order of lim^-jo and 
liniq^o is very important in discussing transport properties since the observable 
currents can be obtained by the dynamical and uniform field (i.e., lim^^o liniq^o) 
but the static and non-uniform field (i.e., liniq_>o lim^^o) does not cause any ob¬ 
servable currents due to the screening induced by the modulations of the charge 
distribution' 3 ^. Also, the value of WTt rans is very important since the adiabatic con¬ 
dition ujT t ra „ s <C 1 means the realization of local equilibrium due to the rapid 
relaxation compared with w -1 , a typical time scale of the field; as a result of that 
relaxation, the electronic transports are governed mainly by the QPs near the Fermi 
level. For example, such importance of the i nequ ality of WTt rans is seen from the dif¬ 
ference between the zero and the first sound' 3 ' 3 . 

Replacing T an d m' 311 ®q. (fl3l) by the corresponding contour 

integrals* 3 ^3 and doing several straightforward calculations^ with attention to 
the analytic properties' 3 ' 3 of the single-particle Green’s function and four-point ver¬ 
tex function, we can carry out the analytic continuation of K a b c d{k,k';ifl n ). As a 
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result, K (0,w) is given by 


E {Vkv)ba{y k ' v )cd 

iV k,k' {a} 

f°° de 
x i 


47T* L tanh 4r K ^ bcd<yk ' fe,; e; ^ 


+ (tanh - tanh K^ bcd (k, k! \ e; u) 


- tanh £ ; w ) 


(14) 


where K iSLd ( k > e ; w ) is 

/ oo / 3 

7 — E E gi\aABb(k\u) 

-00 471-2 j/ = i 

X Jii'-,{A}(k,k';u])gi> ] ccdD{k'-,uj), (15) 

with gi- ac db{k\ui) being 

01 ;a«ft(*;w) = + ( 16 ) 

g2-,aedb(k;u;) = G™(k,e + w)G%\k,e), (17) 

and 

g 3 ;acdb(k-,oj) = G^\k, e + oj)cf\k, e) , (18) 

and Ju'.{A}(k,k';uj) being 


Ju-,{A}(k,k';u>) = tanh— T 11 .i i{A} (k,k';uj) 


+ coth 


e — e 


2 T L 
e' + u 


r 11-II;{A} (k, k'; w) - r n _i ;{ A} (k, k'; u) 


Ji 2 ;{A}(k,k';uj) = (tanh 6 U - tanh r 12 ;{A} (fc, k'\ w), 

^ j 

Ji 3 -,{A}(k,k';u) = - tanh r 13 _ I;{A} (fc, k';u) 


(19) 

( 20 ) 


— coth 


f + e' + wr 


2T 


ri3-ii;{A}(fc,fc';w) - r 13 _i :{j 4}(fc, k'; w)J, (21) 

( 22 ) 


J 2 i-,{A}(k,k';u) = tanh—r 21;{A} (fc,A:';w), 

( 6 ; — e e / \ 

coth - tanh — J r 22 -ii;{ J 4 } (^, w) 


+ ^coth 


2 T 
d + e + uj 
2T 


( tanh L±E_ coth e + e + w 


2 T 


^ j 

^23;{A}(fc,fc';w) = - tanh r 23 ;{A} (fc, fc';a;) 


cot h r 22-III;{A} (k,k’-,U)) 

)r 22 _iv ; {A}(fc,fc / ;w), (23) 

(24) 


2T 
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lm£' lm£— lm£' = 0 


► lm£ 


\ 31-11 

31-1 

21 

11-11 y' 
y' 11-1 


\ 22-11 y' 


32 

22-1 >'22-111 
y'22-IV\ 

12 



‘ s 

33-1 / 

23 

\ 13-1 

/ ' 


y 33-11 


13-11'-. 

lm£ + lmCO = 0 

Im£ + lm 


• lm£'+ lmCO = 0 


Fig. 2 . The connection between the additional subscripts of the four-point vertex function and the 
relations among (e, e',cj) of that function. 


J3i-,{A}(k,k';u) = tanh—r 31 _ I;{A} (/c,fc';w) 


coth ■ 


2 T 


r 3 i-ii;{A}(fc, fe';w) -r 31 _ I;{A} (fc,fc';w)J, (25) 

J32-,{A}(k,k';uj) = (tanh £ + ^ - tanh ^)r 32; { A }(fc, k’\u), (26) 


and 

J33;{A}(k, k'; u) = - tanh ■ 
— coth - 


2 T 


2 T 


- r 33-I ;{A}(M';w) 

r 33 -n -,{A}{k,k'\w) -r 33 _i . {A} (k,k';w) 


(27) 


In Eq. (1151) , I have not explicitly written whether the frequency integral is the prin¬ 
cipal integral or not; the integrals containing hyperbolic cotangent are the principal 
ones. Also, in Eqs. m (GB the additional subscript of the four-point vertex func¬ 
tion such as 12 represents the relations among its three frequency variables, as shown 
in Fig. [2] Since is determined by the Bethe-Salpeter equation, 


JlV; {a} {k,k'-,u) = J l [] ) {a} (k : k';U J )+ fZ^l";a b CD(k,k";U>) 

1 l "=1 iV fc" {A} 


X gi";CABD(k ; u, )'Jl"l’-ABcdJ : >k’,U)), (28) 


and Jii' ) i a i(k, k!\ w) is obtained by the method explained in Sect. 2.3, we can exactly 
calculate a uv from Eqs. m, and (ITU) (12511 in principle. 

Then, to rewrite Eq. (1141) in a more compact form, we use the three-point vector 
vertex function instead of the four-point vertex function. The three-point vertex 
function in Matsubara-frequency representation is defined as 

E GaA(k + q,iui m+n )A V[ AB(k,iu mi Qi 'Bb(fci m) 


A.B 


dre v 


if2 n 


^TCk+qa(j)J~qu{j )£],(,}• 


(29) 


l 


l 
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Since the analytic continuation of the three-point vector vertex function can be 
carried oihE in a similar way to that used for K a i >c( i(k,k , ',i£l n ), we obtain the 
three-point vector vertex function in real-frequency representation, 


Aj s\l\ab(ki Q^) iykv)ab 


f T7 E E E f-oo ztiJll';abCD(k, k'] q) 

iv k' {A}1' = 1 

X gi'-CABD(k'',q)(Vk' v )AB- 


(30) 


Here the additional subscript of the three-point vector vertex function, l, denotes 
the conditions about its e and u>: l = 1 denotes Ime > 0 and line + Imu >0; 1 = 2 
denotes Ime < 0 and Ime + Imw > 0; l = 3 denotes Ime < 0 and Ime + Imcj < 0. 
Combining Eqs. (fTT!) and (fl5l) with Eq. (13U1) . we can rewrite Kf 1 
1 


} (0,u;) 


as 


^)(0,u;)=--EEK,)6a/_ O 


k {a} 


de 
oo 47 xi 


tanh —gi-acdb(k; w)A„ ; i . cd {k; u) 


+ 


^tanh - tanh g 2 -, ac db{k; u)A v - 2 -, c d(k; w) 

- UJ 


- tanh 2r ff 3 ;aedfc(fc; w)A y;3;cd (fc; u) 


(31) 


Since the exact expression is difficult to s olve, we use the most-divergent-term 
approximation^!, based on the properties^ of the product of the single-particle 
Green’s functions in the limit q = (q, to) -A 0 in the presence of the QP peak. This 
is correct in the FL and remains appropriate in a correlated metallic system where 
perturbation expansion is satisfactory. As explained in Sect. 1, applicability of the 
FL theory differs from that of a microscopic perturbation theory. The microscopic 
perturbation theory is satisfactory to describe electronic properties of a correlated 
electron system if perturbation expansion has a good convergence or becomes an 
asymptotic expansion. Let us recall, first, that a well approximate partial sum can 
be constructed even if its convergence is not good^, second, that perturbation 
expansion becomes an asymptotic expansion near a phase transition. 

Befo re p roceeding with the formal derivation of oyj/, we remark on the 
property^ “ of a pair of the single-particle Green’s functions in the limit q —> 0. 
When the QP peak exists and the QP damping is much smaller than T (i.e., cold- 
spot-type), only G^\k + f,e + %)G^\k — f,e — f) gives the most singular 
term being inversely proportional to the QP damping due to the merging of the 
poles of these Green’s functions in q —> 0; the others, the retarded-retarded and 
the advanced-advanced pair, become the higher order terms. In this discussion, we 
have considered only the contribution from the coherent part of the single-particle 
Green’s function since the incoherent part does not lead to such singular term in 
q 0. Also, I have used lim^-^o-i- — a-E-E = — 2.Tti5{x — a) since the QP 

damping at momenta near the Fermi level is assumed to be negligible at low tem¬ 
peratures. The singular contribution from the hot spot is less important than that 
from the cold spot since the QP damping at the hot spot remains non-negligible 
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even on the Fermi level. Since the existence of the QP peak and cold-spot-type 
QP damping is appropriate at least for sev eral momenta in a metallic phase near 
a magnetic QCpE or a Mott transition^, the leading terms in gi ;ac db(k; ui) with 
respect to the QP damping or ui at low temperatures and frequencies are given by 


(7 v—v a \ Za(k)zp(k) 

9l ’ acdb{ ‘ k]U,) [e - £(fc) + *o+][ e _ £*(fc) + *0+] ’ 

a. ^ o_v- n.\ z<*(k)zp(k)6(e - £*(k)) 

92 -^k^) ~ 27T* E u aac , m {k) w _ a(fc) + ^ {k) + f[7 , (fc) + 7 * (fe)] ■ 


and 


g^acdbik-M -Ew(fc) [e _ e * (fc) _ l0+][e _ Cfj{k) _ ,0+] ’ 


(32) 

(33) 


(34) 


where u aac -,dpb{k ) = ( U k ) aa (Ul) ac (U k ) d p(Ul)p b with (U k ) aa being the unitary ma¬ 
trix to obtain the QP bands [not equal to (Uj?) aa ,]. Thus, the most divergent terms 
in the coherent limit arise from g 2 \acdb(k\ uS) at cold spots. 

Using the above property, we obtain an approximate expression of a vv where 
we consider only the most divergent terms with respect to th e Q P damping in the 
coherent limit. To use the most-divergent-term approximation^, we introduce two 
quantities, (k, k'\u) and A ^. ab {k',u>), which are irreducible with respect to 

only a retarded-advanced pair of the single-particle Green’s functions, 


/-r(O) 


(k,k’;u) 


+ E iuE E I-oo ini u U"\abCD 
l "=£2 ^ k" {.4} 


1 


r°° de" rr( 0) 

TO 'JlW-n 


(fc,/c";w) 

x gi"\CABD(k \u)jl„l- A B C d(k ,k-,uj ), (35) 


and 


ASUfc; W ) = (u fc „W+£ E T7E J-ocifiJwUcn^k'^) 

{A} Z'#2 fc' 

X gi'-CABD{k'-,Uj)(vk'iy)AB- 


(36) 


Using these two quantities with Eqs. ®, m- m and (liTOl) and the exchange 
symmetry of the four-point vertex function among its momentum and frequency 
variables, we can rewrite Eq. m as 


kSP M=-^EE M ba fZ £ [tanh ^g 1;acdb (k; u ) A^ ;cd (fc; w) 


k {a} 


47T* - 


k (0) 


2 T" 


- tanh — -g 3 .acdb(k ; w)A^. cd (A;; w) 


( 0 ) 

3;< 


-IfEEA 


IV 


de ( e + uj e \ 

—: I tanh tanh — ' 


k {a} 


4:TTl 


2 T 


2 T) 


x g2-acdb{k;u})A u -2-cd(k- w). 


(37) 


At this stage, this expression remains exact. Then, since only the second term in 
Eq. K37D contains a retarded-advanced pair and the leading term with respect to 
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to comes from (tanh — tanh ss 2u>(— we can obtain an approximate 

expression of <j vl/ in the most-divergent-term approximation, 

de . mi „ „s „ ^ . ,, / df{e)' 


9 3 

= — E E f°° 

AT 1 ' J —c 

h {a}—1 


■A^2;6a( fe ; 0)g 2;acdb (k- : 0)A, ;2;cd (fc; 0) (-^E). (38) 

From this and Eq. (1331) . we can show that oy„ is inversely proportional to the 
unrenormalized QP dampinJSSI. 

A<°> 

J V;2, 

Eq. (1351) with the Ward identitj®, hf\. ba {k] 0) is given by 


2 , ba (k\ 0) and A v - 2t cd(k', 0) i n Eq. (1551) are determined as follows. By combining 


^u;2;ba(^> ®) — { V kv)ba + 


5ReE^(fc) 


V;2;baV ,v J'”V — Iba 1 -• (39) 

In the present model, Re£^(fc) = ReE^(fc) holds due to the even-parity and the 
time-reversal symmetry. Since l , (&,£/; 0) = /o'; 0) is satisfied due to 

Eqs. dm (l26l) . and (l35l) . A u - 2 - c d(k-0) is given by 


A„ ;2;cd (fc; 0) = A^ ;cd (fc; 0) + -EE fZ ^J^cdCD^, V ; 0) 

iv fc' {A} 

x 92 -cabdTT'; 0)A v . } 2;ab(T'-, 0 ). 


1 ^ ^ poo de ' „(1) 


(40) 


Thus, the present theory takes account of the CVCs due to the self-energy 
and irreducible four-point vertex function, neglected in the relaxation-time 
approximation®®. 


2.3. FLEX approximation with the CVCs arising from the 

self-energy and MT term for a multiorbital Hubbard model in 
a PM phase 

To calculate the resistivity, I use the FLEX approximation with the CVCs arising 
from the self-energy and MT term in a PM phase. I first explain the FLEX ap¬ 
proximation for a multiorbital Hubbard model in a PM phase. Next, I derive the 
irreducible four-point vertex function in this approximation; as described in Sect. 
1, I consider only the MT term and neglect the AL term. Then, I show the Bethe- 
Salpeter equation for the current by using the MT term as the kernel of the CVC. I 
also remark on the numerical treatment of the principal integral appearing in that 
CVC. Finally, I discuss applicability of this method. 

I determin e sever al single-particle or two-particle quantities using the FLEX 
approximation^*® SB in a PM phase where only the electron-hole scattering pro¬ 
cesses of the bubble and the ladder diagrams are considered as the Luttinger-Ward 
functiona l®^®** !, $lw[G]. Since that is a conserving approximation®^ ®® based on 
the thermodynamic potential expressed in terms of the single-particle Green’s func¬ 
tion, we can determine single-particle or two-particle quantities by using 3 >lw[G] 
and its functional derivative. Since constructing $lw[G] is equivalent to construct¬ 
ing the effective interaction, where the same kinds of the diagrams are considered, 

















June 17, 2015 0:42 WSPC/INSTRUCTION FILE review-rho'final 


14 Naoya Arakawa 

and the latter is easier, we formulate the FLEX approximation by the latter pro¬ 
cedure as follows. First, considering the electron-hole scattering processes of the 
bubble and the ladder diagrams for H\ nt , we obtain the effective interaction in the 
FLEX approximation, 

1 

2 L~ 

1 

2 L~ 


Vabcd S3Si ( < l' ~ o Wabcd ~ X ^abABXABC Dcd 

{.4} 


Uabed + X ^o6AbXaSCd(9) iftn)Uc 

{A} 


CDcd 


(J® (J® 

^ SlS 2 ^ S 4 S 3 


' S 1 S 2 ^ S 4 .S 3 5 


(41) 


with 


Xabcdi^'^n) — XabcdiQG^n) + X XabAB (Q: ^n)U\ BC dXc Dcd^Q > ^n) > 

Ml 

Xabcd^T^n) = Xabcd{<h i^n) ~ X XabABiQ, ^n)U^ B Q dXc Dcd^ > *^n) i 

{A} 


and 


T 


Xabcd{Q> ifX) — X X ^ac{k 4“ *35 ^m+n)G d b{k^ fkl m ). 


(42) 

(43) 

(44) 


fc m 


Here we introduce the bare four-point vertex functions in spin and charge sector, 


and 


U: 


abed 


= uIL - t/TL, = < 


'abed 


tt _ 


abed 


rrC _ rrtt 
u abcd ~ u - 


abed 


r/tt _ y 

abed * 


u 

for a = b = c = c? 

Jh 

for a = b ^ c = d 

U' 

for a = c b = d 

./' 

for a = ci b = c 

u 

for a = b = c = c? 

to 

q 

for a = b ^ c = d 

-t/' + 2 J H 

for a — c ^ b = d, 

J' 

for a = d ^ b = c 


(45) 


(46) 


where t/ 0 s b s cd is t/ a s | cd = £/ a s b s c s / and the spin-flipping term, C/ 0 s b s c f for s ^ s', sat¬ 
isfies = — t^dcM = — ^abed * Also, we neglect the vertex corrections to the 

susceptibilities in spin and charge sector; its effects are discussed later. Then, the 
single-particle Green’s function is determined by the Dyson equation, 

G a 5(/c, ZCC m ) — G ab {k, iuJm) 4“ G a A (^, i'Ldm)^- 1 AB (^, i^rn)G Bb (^, i^m) , (47) 

A.B 


with the self-energy given by 

T 


^ac{k^ iuJm) — X X X VabcdiQi 'l n )G b d{k Qi ^m—n)i 

t* q n b,d 


(48) 
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where 


Vabcd{q,i£ln) — iO, n ) V^bcdfai ^n) S UaAbBXABC D(<h i®"n)U^f cDd 

{A} 

3 r i 

Uabcd d(Q' i^nWcDall 


2 L 




1 

+ 2L 


Uabcd + J2 ^abABXABCD(.q^^ri)U§ Dcd 
{. 4 } 


UaAbBXABCD(q, iftn)U(f cDd . (49) 

MI 


The last term in Eq. dm is introduced to exclude the double counting of the 
topologically equivalent term in the self-energy. Solving Eqs. (El) (Ell and (S3 
(El) with Eqs. ((HI)' (1TU1) . (H5l) and (El) selfconsistently by iteration, we can determine 
the single-particle and the two-particle quantities in the FLEX approximation. 

It should be noted that the partial inclusion of mode-mode couplings^ for 
fluctuations, the interactions between fluctuations at different momenta, through 
the self-energy improves some unrealistic results obtained in the random-phase 
approximation, although the susceptibilities are determined by the random-phase 
approximatio n-type (but renormalized) equations. For example, in the FLEX 
approximation^® 2 9|61 1 65 1 va j ue 0 f jj f or a magnetic transition becomes about 2 

eV, the momentum dependences of the mass enhancement and FS deformation are 
taken into account, and the Curie-Weiss-like temperature dependence of the spin 
susceptibility is obtained near a magnetic QCP. In particular, the final improvement 
is powerful to describe electronic properties near a magnetic QCP. 

We also determine the irreducible four-point vertex function in the FLEX ap¬ 
proximation in keeping conservation laws^U In a conse rving approximation, the 
irreducible four-point vertex function is given bjEH 31 63j 


■p(l) 

abed 


(fc, 


k',iu m r,q,iO, n ) = 


ST, ab (k,iu} 

m) 

SGcd(k', iu) m ') ’ 


(50) 


For the actual calculations, we first calculate the right-hand side at q = 0 and 
Q n = 0 and then label momentum and frequency transfers correctly as the electron- 
hole scattering process among an electron of orbital b with (k,iui m ), a hole of 
orbital d with (k 'an electron of orbital a with (k + q,iui m +n), and a hole 
of orbital c with ( k' + q,iuj m >+ n ). After several straightforward calculationsPl by 
using Eqs. (l48l) (f50jl . we obtain the irreducible four-point vertex function in the 
FLEX approximation, which is the sum of the MT and the AL ternPSl. In this 
paper, I consider only the MT term, 


AL 

' abed 


(fc, iid-m-, & 5 tUJ-m' \ Q-, — Vacbd{}z ^ 


(51) 


This treatment will be sufficient for a qualitative discussion about many-body effects 
on p a b since the CVC arising from the AL ter m g ives the higher order contribution 
compared with that arising from the MT termed I have checked the validity of this 
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statement by calculating the main terms of the AL CVC . Since it is necessary to 
calculate J-nl-dCD^- ^ 0) in Eq. (l40l) . we need to carry out the analytic continu¬ 
ation of Eq. (1FT1) in region 22-11, region 22-III, and region 22-IV. Carrying out the 
analytic continuation^ of the MT term and using Eq. (l23l) . we obtain 

*^22 }cd,CD ( k , k'- 0) = 2i(coth < -^ r + tanh ^)lmV c ( *] D (fc - k'). (52) 

Substituting Eq. (15211 into Eq. (FTTT1) . we obtain the Bethe-Salpeter equation for 
the current in the FLEX approximation with the MT CVC in a PM phase, 


A v,2;cd(k; 0) = A i% iCd (k) + ^EE/“| ( coth e ~YF + tanh i) 

k' {A} 

x Im V^\ D (k - k')g 2 -cABD{k 0)A y;2 -AB{k'\ 0). (53) 

We see the roles of the MT CVC are similar to those of the backflow correction^. 

Before discussing applicability, I explain how to treat the principal integral in 
Eq. (1531) for the numerical calculations. Since both the numerator and denominator 
of the term containing coth in the MT CVC become zero simultaneously at 
e' = e due to Im V^ b J d (q, 0) = 0, the principal integral can be calculated as follows: 


/ OO 7 / / 

^ coth ^-^-Im V cCd D ( k - k')g 2 -cABD(k 0)A v - 2 . A B(k; 0) 

= / tt - coth ^-Im V^d D (k - k')g 2 - C ABD{k' ; 0)A v - 2; AB(k; 0) 
^^T r ^7[( e ^ + l ) lmV cCdDi k - k')g 2 -cABD{k']Q)K v . 2 .AB{k'-,Q)\ 


(54) 


where the hrst term contains the contributions other than e' = e. 

Finally, I discuss applicability of the FLEX approximation with the CVCs aris¬ 
ing from the self-energy and MT term. 

First, the FLEX approximation is suitable to describe the electronic structure 
at low temperatures for moderately strong electron correlation. In a single-orbital 
Hubbard model on a square latticd^22l, the imaginary-time dependence of the 
single-particle Green’s function at several momenta in the FLEX approximation 
shows satisfactory (but not perfect) agreement with that in the quantum-Monte- 
Carlo calculation at U being a half of the bandwidth; the agreement becomes better 
near the AF QCP than away from the AF QCP. Since the similar agreement will 
hold even in a multiorbital Hubbard model on the same lattice and the FLEX ap¬ 
proximation can treat the coherent part of the electronic spectrum satisfactorily ^21, 
the electronic structure in metallic phases of the present model at low temper¬ 
atures will be well described by the FLEX approximation at least qualitatively. 
Actually, the FLEX approximation succeeded in reproducing the larger effective 
masJ^S of the d xy orbital than that of the d xz / yz orbital^® and its agreement with 
the experiment®®] is better than the case of the DMFT^l, as described in Sect. 1. 
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In contrast, the FLEX approximation becomes unsuitable at high temperatures 
or near a Mott transition for strong electron correlation. This is because local 
correlation plays important roles in such casi^® ^® and the effects of local correlation 
on the electronic spectrum are smeared out in the FLEX approximation^®. 

Then, the magnetic properties at low temperatures for moderately strong elec¬ 
tron correlation will be appropriately described by the FLEX approximation if the 
dominant correlation of the system is spin fluctuation whose largest contribution 
comes from a non-degenerate orbital. Due to neglecting the vertex corrections to 
the susceptibilities, the enhancement of spin fluctuation arising from electron cor¬ 
relation is overestimated in the FLEX approximation compared with the enhance¬ 
ment of charge or orbital fluctuation. Actually, in a two-degenerate-orbital Hub¬ 
bard model on a square lattice at small (Jh//7) near an AF QCP, the AL vertex 
correction to the susceptibilities causes the enhancement of orbital fluctuation®®. 
However, I believe that in the present model the FLEX approximation is sufficient 


ualitatively since the d xy orbital gives 
in this case, even if the MT and 


to describe the magnetic properties at least 
the largest contribution to spin fluctuatio: 
the AL vertex correction to the susceptibilities are considered beyond the FLEX 
approximation, the magnetic properties will not qualitatively change and these 
corrections will modify the values of the susceptibilities since orbital fluctuation 
enhanced due to the AL term does not dominate over spin fluctuation. Actually, 
the strongest enhancement of spin fluctuation at Qic-af = (^r> Tjf 1 ) i n the FLEX 
approximation^ away from and near the AF QCP agrees with the experiments in 
S^RuC)^® and S^Ruo. 075Tio,025 oJ3, respectively. 

Moreover, if the vertex corrections to the current arising from the self-energy 
and irreducible four-point vertex function due to electron correlation are added to 
the FLEX approximation, this method is suitable to describe the transport prop¬ 
erties of a metallic phase due to low-frequency external field satisfying W7t ra ns "C 1 
at low temperatures. In contrast to case of the vertex corrections to the susceptibil¬ 
ities, the vertex corrections to the current are essential for discussing the transport 
properties^ since the CVCs are vital to satisfy conservation laws and conservation 
laws play significant roles in transport phenomena. For example, the importance of 
the treatment holding conservation laws is known for a system without the lattice 
(e.g., the electron gas): only if the CVCs due to electron correlation are correctly 
taken into ac cou nt, we can obtain the correct results such as the absence both of 
the resistivity®® and of the renormalization of the Drude weight®'' 7 and electron cy¬ 
clotron frequency®®. Also, the CVCs due to electron correlation are important in a 
system with the lattice since these CVCs are necessary to obtain the correct effects 
of electro n corr elation on the transport coefficients in the presence of the Umklapp 
scattering®®®'' 7 ®®. Another example showing the importance of the CVCs is the 
emergence of the Curie-Weiss-like temperature dependence of the Hall coefficient 
near the AF QCP due to the MT CVC in the FLEX approximation in the single¬ 
orbital Hubbard model on a square lattice^®. In addition to the treatment holding 
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(a) (b) (c) 



Fig. 3. Momentum dependences of the static spin susceptibility, x s (q, 0) = YL a h m,(9> 0), 
for several temperatures at U = (a) 1.8 and (b) 2.1 eV, and (c) temperature dependences of 
X S (QlC-AF,0) at U = 1.8 and 2.1 eV. 


conservation laws, the satisfactory treatment of the coherent part in the FLEX 
approximation is powerful to describe the transport properties of a metallic phase 
due to the low-frequency external field. This is because the dominant contributions 
to the response induced by that external held at low temperatures come from the 
contributions near the Fermi level as a result of the energy derivative of the Fermi 
distribution function in the response function. Furthermore, this powerfulness of 
the FLEX approximation will hold even near a Mott transition since it is shown in 
the DMFT-21 for a single-orbital Hubbard model that the transport properties in 
wTt rans -C 1 are well described in the approximation where only the coherent part 
is considered. Thus, the FLEX approximation with the CVCs arising from the self¬ 
energy and MT term is satisfactory to describe the transport properties of metallic 
phases of the present model at low temperatures. Actually, this method^ near and 
away from the AF QCP reproduced the temperature dependence of several trans¬ 
port properties of Sr 2 Ru 04 ^ and Sr 2 Ruo.o 75 Tio. 02504 ^, as pointed out in Sect. 
1 . 


3. Results 

In this section, I show the results of p a b = &xx °f some quasi-2D PM ruthenates in 
the FLEX approximation with the CVCs arising from the self-energy and MT term 
and compare these results obtained in more simplified cases than that method. In 
particular, we focus on the effects of the self-energy and MT term due to electron 
correlation and the role of each t 2ff orbital. 

I carried out the numerical calculations as follows. I set the 64 x 64 meshes of 
the Brillouin zone and 2048 Matsubara frequencies and used the fast Fourier trans¬ 
formation with the zero padding method . I obtained the single-particle Green’s 
function, self-energy, and MT term in the FLEX approximation by solving Eqs. 
11321) (1331) and (1371) (1351) with Eqs. ©, (133)1) . (1351) and (1351) by iteration, where I as¬ 
sumed that convergence was reached when the difference between the self-energies 
before and after certain iteration was less than 10 -4 . To obtain the quantities as 
a function of real frequency, I used the Pade approximation^ using the quanti¬ 
ties at the lowest four Matsubara frequencies. The real-frequency integrations were 














June 17, 2015 0:42 WSPC/INSTRUCTION FILE review-rho'final 


Many-body effects on the resistivity of a multiorbital system beyond Landau’s Fermi-liquid theory 19 



T(eV) 


Fig. 4. Temperature dependence of p a b at U = 1.8 eV. The definition of each case is described in 
the main text. The inset shows p a b against T 2 below T = 0.01 eV. 


approximated by the integrations with the interval 0.0025 eV and the upper and 
lower values 0.2 and —0.2 eV. The current was determined by solving Eq. (l53l) by 
iteration, where the convergence condition was assumed to be that the difference 
between the currents before and after certain iteration was less than 10~ 4 . 

In the following, I consider cases at U = 1.8 and 2.1 eV as cases away from and 
near the AF QCP, respectively. These correspondences are because of the Pauli-PM 
temperature dependence of the spin susceptibility at U = 1.8 eV and the Curie- 
Weiss-like temperature dependence of the spin susceptibility at q = Qic-af at 
U = 2.1 eV [see Figs. GD[a)-[3][c)]. The latter is characteristic of a magnetic QCP 
and causes the strong temperature dependence of the Landau parameters with 
momentum transfer Qic-af through the temperature dependence of the reducible 
four-point vertex function. The choice of U i s reasonable since the value estimated 
experimentally in Sr 2 Ru 04 is about 2 e\m 

As explained below in detail, there are three main results: (i) p a b of some quasi- 
2D PM ruthenates without the rotation and the tilting of RuC>6 octahedra is de¬ 
termined almost by the conductions of the d xz / yz orbital due to the smaller un¬ 
renormalized QP dampings of the d xz / yz orbital than those of the d xy orbital; (ii) 
The crossover between the T-linear and the T 2 dependence of p a b occurs away from 
the AF QCP at about T = 0.008 eV due to the temperature dependences of the 
unrenormalized QP dampings of the d xz / yz orbital; (iii) The T-linear p a b emerges 
near the AF QCP due to the hot-spot structure of the QP dampings of the d xz / yz 
orbital at k = Qic-af and its symmetrically equivalent ones. 

3.1. Case away from the AF QCP 

I begin with the temperature dependence of p a b at U = 1.8 eV in Fig. Q] Three 
cases in that figure are defined as follows: in MT CVC case, the CVCs arising from 
the self-energy and MT term are included; in No MT CVC case, the CVC arising 
from the self-energy is included; in No all CVCs case, which is equivalent to the 
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(a) (b) (c) 



Fig. 5. Temperature dependences of the total and orbital-decomposed components of cr xx at 
U = 1.8 eV in (a) the MT CVC case, (b) the No MT CVC case, and (c) the No all CVCs case. 
The definition of each orbital-decomposed cr xx is described in the main text. 


(a) (b) (c) 



Fig. 6. Momentum dependence of (a) the unrenormalized QP damping and (b) the QP damping for 
each t 2 g orbital at (U, T) = (1.8, 0.01) (eV), and (c) temperature dependence of the unrenormalized 
QP damping of the d xz / yz orbital at k = Q ic-AF at U = 1.8 eV. The dashed line in panel (b) 
denotes T = 0.01 eV. The inset in panel (c) shows the data against T 2 . 


relaxation-time approximation^, all the CVCs are neglected. There are three main 
remarks about Fig. 2) First, in all the three cases, decreasing temperature causes the 
crossover from the T-linear dependence to the T 2 dependence at about T = 0.008 
eV. The T 2 dependence at low temperatures is more clearly seen from the inset 
of Fig. U Second, the value of p a b in the No all CVCs case at each temperature 
is largest in the three cases. Third, the value of p a b increases from that in the No 
MT CVC case when the MT CVC is included. The first remark indicates that the 
CVCs little affect on the power of the temperature dependence of the resistivity. The 
second indicates that the value of the resistivity is overestimated in the relaxation¬ 
time approximation. The third indicates that the MT CVC enhances the resistivity 
as a result of the reduction of the current, which is similar to the effect of the 
backflow correction. Thus, the main effects of the CVCs on the resistivity just 
change its value. 

Then, to understand the role of each t 2 g orbital, I analyze the orbital- 
decomposed <j xx at U = 1.8 eV. The orbital-decomposed a xx for the d xz / yz and the 
d xy orbital are obtained by replacing ^, o > =1 in Eq. (1551) by ]C{a}=i and Xq a }= 3 > 
respectively. Those components are sufficient in the present model since the intraor¬ 
bital components are much larger than the interobital components due to the large 
intraorbital hopping integrals compared with the interorbital ones. We see from 
Figs. EJa)-[5][c) that the dominant contributions to the total of <j xx come from the 
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T(eV) 


Fig. 7. Temperature dependence of p a b at U = 2.1 eV. The definition of each case is described in 
the main text in Sect. 3.1. 

component of the d xz orbital, and that the component of the d xy orbital is less than 
10% of the total. Due to the rotational symmetry of the system, in case of (J yy , the 
component of the d yz orbital gives the dominant contributions. Thus, the inplane 
transport is governed mainly by the conductions of the d xz / yz orbital. 

The above orbital-dependent transport arises from the smaller unrenormalized 
QP dampings of the d xz / yz orbital than those of the d xy orbital since ay v is inversely 
proportional to the unrenornralized QP damping, as explained in Sect. 2.2. Actually, 
we see from Fig. [HJa) that the unrenormalized QP dampings of the d xz / yz orbital 
are smaller. The similar orbital dependence holds at the other temperatures. 

Also, we see from Figs. Ob) and[6])c), respectively, that the QP dampings of the 
dxz/yz orbital remain the cold spot, and that the temperature dependence of the 
unrenormalized QP damping of the d xz / yz orbital at k = Qic-af changes from T 2 
to T-linear at about T = 0.008 eV. Thus, the latter is the origin of the crossover of 
the power of the temperature dependence of p a b at about T = 0.008 eV. 

3.2. Case near the AF QCP 

I turn to the temperature dependence of p a b at U = 2.1 eV in the three cases, 
considered in Sect. 3.1. From Fig. [71 we see that the T-linear p a b emerges in all the 
three cases. We also see the similar effects of the CVCs on the value of p a b to those 
at U = 1.8 eV, the overestimation in the relaxation-time approximation (i.e., the 
No all CVCs case) and the increase from the value in the No MT CVC case due 
to the MT CVC. The first result indicates that the emergence of the T-linear p a b 
near the AF QCP arises from the temperature dependence of the unrenormalized 
QP damping. Furthermore, those results and the corresponding results at U = 1.8 
eV suggest that the effects of the self-energy and MT term on the value of p a b and 
power of the temperature dependence of p a b are ubiquitous. 

Also, comparing Fig. [7] with Fig. [TJ we see electron correlation enhances the 
value of p a b at each temperature. This is due to an increase in the unrenormalized 
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(a) (b) (c) 



Fig. 8. Momentum dependences of (a) the unrenormalized QP damping and (b) the QP damping 
for each t 2 g orbital at ( U,T ) = (2.1,0.01) (eV), and (c) the temperature dependence of the 
unrenormalized QP damping of the d xz / yz orbital at k = Q ic-AF against T 0,5 at U = 2.1 eV. 
The dashed lines in panel (b) denotes T = 0.01 eV. 


(a) (b) (c) 



Fig. 9. Temperature dependences of the total and orbital-decomposed components of (T X x at 
U = 2.1 eV in (a) the MT CVC case, (b) the No MT CVC case, and (c) the No all CVCs case. 
The definition of each orbital-decomposed a xx is described in the main text in Sect. 3.1. 


QP dampings of the d xz / yz orbital with increasing U [e.g., see Figs.^a) and|5Ja)]. 

Then, in a similar way to that used at U = 1.8 eV, I analyze the role of each t 2 g 
orbital at U = 2.1 eV. Figures[5Ja) ]5[c) show the temperature dependences of the 
total and orbital-decomposed components of a xx at U = 2.1 eV in the three cases. 
The orbital-dependent transport holds even near the AF QCP: the conductions 
of the d xz / yz orbital mainly contribute to the inplane transport. Its mechanism is 
the same as that at U = 1.8 eV, i.e. the smaller unrenormalized QP dampings of 
the d xz / yz orbital than those of the d xy orbital [see Fig. [Ha)]. Thus, this orbital- 
dependent transport is characteristic in some quasi-2D PM ruthenates. 

Moreover, we find from Figs. [Hfb) and[8][c), respectively, that the QP damping 
of the d xz / yz orbital at k = Q ic-af becomes the hot spot, and that the temperature 
dependence of the unrenormalized QP damping of the d xz / yz orbital at k = Qic-af 
is roughly proportional to T 0 - 5 . Since such T 0 - 5 dependence near the hot spot causes 
the T-linear depen den ce of the average of the unrenormalized QP dampings for 
states along the Fl£E2, the origin of the T-linear p a b near the AF QCP is the hot¬ 
spot structure of the QP dampings of the d xz / yz orbital at momenta connected each 
other by the characteristic spin fluctuations of this AF QCP. 
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4. Conclusions 

In summary, I reviewed many-body effectd^S of p a b of the quasi-2D PM ruthenates 
away from and near the AF QCP in the FLEX approximation with the CVCs arising 
from the self-energy and MT term or with the CVC arising from the self-energy or 
without all the CVCs. 

The temperature dependence of p a b away from and near the AF QCP quali¬ 
tatively agree with experiments of Sr 2 Ru 04 ' and Sr 2 Ruo.o 75 Tio .025 04^1 respec¬ 
tively: in case away from the AF QCP, the crossover between the T-linear and the 
T 2 dependence at about T = 0.008 eV and the T 2 dependence at low temper¬ 
atures are obtained; in case near the AF QCP, the T-linear dependence even at 
low temperatures is obtained. Here the main effect of the Ti substitution on the 
temperature dependence of p a b is assumed to be pushing the system nearer the AF 
QCP than Sr 2 Ru 04 . 

The obtained results reveal some important aspects of many-body effects on the 
resistivity of correlated electron systems. First, the overestimation of the value of 
Pab in the relaxation-time approximation and the back-fiow-like effect of the MT 
CVC on the value of p a b are ubiquitous. It is also ubiquitous that the power of the 
temperature dependence of the resistivity is determined almost by the temperature 
dependence of the momentum- and orbital-dependent unrenormalized QP damping. 
The similar results have been obtained in a single-orbital Hubbard model on a 
square lattioPHI Moreover, the T-linear resistiv ity near the AF QCP is similar 
to that obtained in those previous studies^ ^ However, I emphasize that the 
criticality of the resistivity, the power of its temperature dependence, does not 
always connect with the criticality of spin fluctuation enhanced near a magnetic 
QCP in multiorbital systems, while these are always the same in single-orbital 
systems. This characteristic property comes from the facts that the orbital whose 
unrenormalized QP damping is small mainly contributes to the resistivity, and that 
the main orbital of the characteristic spin fluctuation of a magnetic QCP has the 
large unrenormalized QP damping. Since the momentum, the temperature, and the 
orbital dependence of the damping of a QP (i.e., the unrenornralized QP damping 
or the QP damping) are overlooked in Land au’s FL theory and such momentum 
dependence is overlooked in the DMFT 1511SI the obtained results highlight their 
importance in discussing the resistivity of correlated electron systems. 

I close this paper with several remarks about the remaining issues . First , it is 
necessary to study the transport properties of other ruthenateJ^ — — using 
the method I used and discuss the similarities and differences. In particular, the 
study for a quasi-2D ruthenatd^ near a ferromagnetic QCP is highly desirable to 
understand the similarities and differences between many-body effects and role of 
each f 2g orbital near the AF and the ferromagnetic QCP. Furthermore, it is im¬ 
portant to analyze the transport properties of 3D ruthenates^ since comparison 
of the results in quasi-2D and 3D ruthenates leads to a deep understanding of the 
dimensionality. Other remaining issues are the applications to other correlated elec- 
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tron systems such as transition-metal oxidetfSI, organic conductors^, CeCoInJ^, 
and UPt 3 79 since these studies are important to deduce ubiquitous properties of 
correlated electron systems and characteristic properties of multiorbital systems. 
Then, it is intriguing to study the transport properties in a superconducting phase 
by extending the present method^ in a PM phase since in some caseiP^l the CVCs 
arising from not only spin fluctuations but also superconducting fluctuations play 
important roles. Furthermore, another remaining issue is to clarify the role of each 
t 2 g orbital in the thermal transport® 9 in the superconducting phase of Sr 2 Ru 04 on 
the basis of the method where the orbital dependence of the dampin g of a QP is 
satisfactorily consid ered. T his is because the combination of my result 9 and sev¬ 
eral previous studies®^ ®^ suggest the existence of the difference between the main 
orbitals of the inplane transport and the superconductivity; if this is correct, we 
should pay attention to the effects of the orbital-dependent damping of a QP on 
the thermal transport® 9 -! for correct understanding of the experimental results. 
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